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Abstract 

We extend the dipole formalism for massless and massive partons to random polarisations 
of the external partons. The dipole formalism was originally formulated for spin-summed 
matrix elements and later extended to individual helicity eigenstates. For efficiency reasons 
one wants to replace the spin sum by a smooth integration over additional variables. This re- 
quires the extension of the dipole formalism to random polarisations. In this paper we derive 
the modified subtraction terms. We only modify the real subtraction terms, the integrated 
subtraction terms do not require any modifications. 



1 Introduction 



The recent years have witnessed significant progress in the field of next-to-leading order (NLO) 
calculations. In particular through new techniques for the computation of the virtual one-loop 
corrections a description of multi-parton final states at NLO is no longer out of reach. Multi-jet 
final states are a typical signature in the current LHC experiments and a calculation at NLO is 
mandatory for an accurate and precise description of these signatures. 

The improvements of the techniques for the computation of the virtual one-loop corrections 
occured along three lines: First of all, the traditional approach based on Feynman graphs has 
been brought to perfection [fTT-fTTI. Secondly, techniques based on unitarity or cuts have made a 
major contribution to the field [[T4I - I3TI . Finally purely numerical methods for the computation 
of the virtual one-loop part were shown to be competitive in the multi-parton domain Il32l - l38l . 
With the help of these techniques several multi-parton processes have been computed [|39] - [52l . 

In view of these improvements for the computation of the virtual one-loop corrections one 
has to reconsider potential bottle-necks which hinder us from going to even higher parton multi- 
plicities. In this re-evaluation the real contribution comes back into focus. For a NLO calculation 
the real contribution is given by the Born amplitude squared with one additional parton. The in- 
tegration over the phase space of one unresolved parton is divergent. To render this integration 
finite, the subtraction method is usually employed Il53l - l68l . In current calculations it turns out 
that the subtracted real contribution is the most challenging one, not from a conceptional point 
of view, but in terms of CPU time. 

Modern calculations are often based on helicity amplitudes. The use of helicity amplitudes 
avoids the A^ iagram -scaling of the squared Born amplitude, which would occur if the Born am- 
plitude is expressed in terms of Anagram Feynman diagrams. Instead one sums over all helicity 
configurations explicitly. By using explicit representations for the polarisation vectors and the 
polarisation spinors the amplitude can be calculated as a complex number for any given helic- 
ity configuration. Taking the norm of a complex number is a simple operation and negligible in 
terms of CPU time. This reduces the complexity from A^ iagram down to A/diagram- Furthermore, the 
Anagram-scaling can be avoided by calculating the amplitude through recurrence relations [|69l . 
Let us now discuss the costs of the helicity method: Both quarks and gluons have two spin states. 
It follows that for an amplitude with n external partons we have a sum over 2" helicity configu- 
rations. It is well known that for the Born contribution and the virtual contribution the sum over 
2" helicity configurations can be replaced by an integration over n random polarisations. This 
integration can be done by Monte Carlo methods and can be combined with the Monte Carlo 
integration over the phase space of the final state particles . Doing so will speed up the evaluation 
of the integrand by a factor of up to 2". (The speed up can be less than 2", if intermediate results 
from one helicity configuration are cached and re-used for another helicity configuration.) As 
a slight disadvantage we have to take into account that we integrate now over more dimensions 
(phase space plus n helicity variables) and that the integrand fluctuates more strongly. This im- 
plies that a larger number of integrand evaluations is needed to reach a prescribed Monte Carlo 
integration error. Despite this slight draw -back it turns out empirically that there is a net gain by 
using the Monte Carlo integration over the random polarisations as compared to summing over 
all helicity configurations. 
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It is therefore natural to extend this method to the real contribution. Here, a technical prob- 
lem has to be circumvented: As mentioned above we need subtraction terms in order to render 
the integration over the phase space of one unresolved particle finite. The subtraction terms are 
constructed to match locally the singular behaviour of the integrand. Currently known subtrac- 
tion terms, like the ones within the dipole formalism [|55l - l59l . match locally the behaviour of the 
spin-summed amplitude squared. However, they do not match locally the singular behaviour of 
a squared amplitude with random polarisations. In this case, additional subtractions terms are 
needed. In this paper we provide these additional subtraction terms. The additional subtraction 
terms, which we derive, have the following properties: 

1 . The additional subtraction terms integrate to zero. This implies that we only need to modify 
the unintegrated subtraction terms, which are subtracted from the Born amplitude squared 
with (n+ 1) partons. The integrated form of the subtraction terms, which is added back, is 
not modified. 

2. The additional subtraction terms are independent of the subtraction scheme which is used 
for the spin-summed squared amplitude. They can be used with any subtraction scheme, 
which provides local subtraction terms for the spin-summed squared amplitude. In this 
context we remark that the currently known antenna subtraction terms do not provide local 
subtraction terms for the spin-summed squared amplitude. The antenna subtraction terms 
approximate only the angular average of spin-summed squared amplitudes and are thus not 
local. Antenna subtraction has to be used with an additional slicing parameter. 

3. The additional subtraction terms are formulated in such a way, that they are independent 
of the actual definition of the polarisation vectors and polarisation spinors. We give a 
representation for the polarisation vectors and polarisation spinors in the appendix, but the 
reader may use his own definition. 

We would like to point out that recently local subtraction terms for squared amplitudes with 
helicity eigenstates have been derived 117011711 . These subtraction terms are however not sufficient 
for the method of random polarisations. This point will be discussed in detail in section[2] 

This paper is organised as follows: In the next section we introduce the framework. We 
review the subtraction method and the method of random polarisations. In section [3] we discuss 
the soft and collinear limits in the case of random polarisations. Section 0] contains the main 
results of this paper and gives the additional subtraction terms needed to use the method of 
random polarisations for the real contribution. Finally, section \5\ contains our conclusions. In an 
appendix we give a choice for the polarisation vectors and polarisation spinors. 
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2 The framework 



The starting point for the calculation of an infrared safe observable O in hadron-hadron collisions 
is the following formula: 

(O) = £ / dxif a (xi) [ dx2f b {x 1 )—— — -Jr j- — (1) 

TLJ J 2K(S K Din 1 Kmn 2 ^colour 1 "colour (2 



ab J J " spin VV 'spin 

J d$ n (pi,p 2 ;p3,--, Pn+2 )0(pi,..., p n+2 ) | A,i+2 \ 2 ■ 



(2) 



This equation gives the contribution from the n-parton final state. The two incoming particles are 
labelled p\ and p 2 , while p^ to p n+ 2 denote the final state particles. f a (x) gives the probability of 
finding a parton a with momentum fraction x inside the parent hadron h. 2K(s) is the flux factor, 
for massless partons it is given by 2K(s) = 2s. The quantity n S pin(0 denotes the number of 
spin degrees of freedom of the parton i and equals two for quarks and gluons. Correspondingly, 
"colour (0 denotes the number of colour degrees of freedom of the parton i. For quarks, this 
number equals three, while for gluons we have eight colour degrees of freedom. Dividing by the 
appropriate number of degrees of freedom in the initial state corresponds to an averaging. d§ n is 
the phase space measure for n final state particles, including (if appropriate) the identical particle 
factors. The matrix element |-^+2| 2 is calculated perturbatively. 

The contributions at leading and next-to-leading order are written as 

(0) L0 = Jo n do B , 

n 

(0) NL0 = I O n+l do R + Jo n do v + Jo n do c . (3) 

n+l n n 

Here a rather condensed notation is used. do B denotes the Born contribution, whose matrix 
elements are given by the square of the Born amplitudes with (n + 2) partons l-^+^l 2 - Similarly, 
do R denotes the real emission contribution, whose matrix elements are given by the square of the 
Born amplitudes with (n + 3) partons |-^+ 3 1 2 . do v gives the virtual contribution, whose matrix 
element is given by the interference term of the one-loop amplitude with (n + 2) partons 

with the corresponding Born amplitude ^5 2 - denotes a collinear subtraction term, which 
subtracts the initial-state collinear singularities. Taken separately, the individual contributions 
are divergent and only their sum is finite. In order to render the individual contributions finite, 
such that the phase space integrations can be performed by Monte Carlo methods, one adds and 
subtracts a suitably chosen piece [|55l - l59l : 

(0) NL0 = J (o n+l do R -O n da A ^+jlo n do v + O n do c + O n jdo A \. (4) 

n+l n \ 1 / 

The squared matrix elements involve a sum over the spins of the external particles. Since both 
quarks and gluons have two independent spin states, we can label the independent spin states of 
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QCD partons by "+" and "— ". We thus have for a tree-level amplitude with n external partons 



^ 0) (pi,...,Pn) 



A,] ,...,A,„ 



(5) 



with %i G {+,—}. The amplitudes (pi,Xi,...,p n ,X n ) are called helicity amplitudes. The 

spin-summed squared matrix element l-^ ^ 2 involves therefore the evaluation of 2 n helicity 
amplitudes. For multi-parton amplitudes it is desirable to avoid the exponential growth due to 
the spin sum. One possibility is to replace the spin-summed squared matrix element by a random 
sampling of squared helicity amplitudes. In combination with the phase space integral, which 
usually is done by Monte Carlo techniques, one would therefore write 

2 



^ 2 (p U ..., Pn+2 ) 



(6) 



J d n+2 U j d§ n O n ^%{piM u l)i---iPn+2,Hu n +2)) 



ro,ii«+ 2 



with 



X{u) 



+ for 0<m<4, 



for 



< u < 1. 



(7) 



We refer to this method as "helicity sampling". Within this method, the integrand as a function 
of the variables u\, u n+ 2 is discontinuous at ui = 1/2. 

A second possibility avoids this discontinuity and is referred to as the method of "random 
polarisations". We discuss this method first for gluons: Instead of polarisation vectors £+ and £^ 
with definite helicity "+" and "— " one introduces linear combinations [1721 



e„(4>) = e'^ + e-%. (8) 
We may then replace the summation over the spin states by an integration over the angle (]). 



2k 



x=± 



(9) 



We thus arrive at 
d$ n O n 



■.Pn+2, 



d" 2 II J d§ n O n fl^ 2 (pi,§i,...,p n+2 ,tyn+2) 



[0,11"+ 2 



(10) 



with §i = 2%m. The integrand is now a smooth function of the additional variables The 
method is easily extended to quarks. One defines the linear combinations 



v(4>) 



(11) 
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Then 



271 271 

E u X u X = ^Jd* u(fau(fa, £ v¥ = ±Jd$ v((b)v((b). (12) 

A,=± o X=± 

We have outlined the method of random polarisations for the Born contribution. The method of 
random polarisations works without modification for the virtual contribution and the insertion 
term as well. However, it will not work without modifications for the subtracted real contri- 
bution. The reason is the following: The dipole subtraction terms are constructed by requiring 
to match locally the singular behaviour of the spin-summed squared matrix element of the real 
contribution. If we consider instead individual random polarisations, there is no guarantee that 
the subtraction terms match locally the singular behaviour, in general they will do so only after 
integration over the additional angles. Therefore, if one wants to use the method of random po- 
larisations for the subtracted real contribution, one has to add additional subtraction terms which 
ensure that the integrand is locally integrable. We note that in ref. [17011711 subtraction terms 
for helicity eigenstates have been derived. These subtraction terms allow to use the method of 
helicity sampling for the real contribution, but not the method of random polarisations. This can 
be seen as follows: For example, for the gluon polarisation vectors one has 



e^)*evO) 



e-%+* + e%~ 



Ref. irTmiTTTl provides local subtraction terms for the contributions proportional to £+*£^ and 
£ p *£y, but not for the ones proportional to the mixed combinations £^*£y~ and £^ *£+ . In this 
paper we derive the local subtraction terms required to combine the dipole subtraction scheme 
with the method of random polarisations. We note that in eq. (fT3l) the terms 



% e v (14) 

on the right-hand side are just the usual spin sum. These are approximated by the spin-summed 
dipole subtraction terms. In combining the dipole subtraction scheme with the method of random 
polarisations we need in addition subtraction terms corresponding to the terms 

e- 2 'h+*£- and e 2i %*£+. (15) 

In an individual dipole subtraction term corresponding to the splitting (ij) — >■ i + j with spectator 
k we will have to include the dependence on the helicity angles fa and fyj of the particles i and 
j. The helicity angle fa of the spectator does not enter the dipole subtraction terms. The terms 
independent of fa and §j will just be the spin-summed dipole subtraction terms, which are well- 
known and given in ref. Il5~5ll56ll59l . Therefore we only have to derive the dipole subtraction 
terms with a non-trivial dependence on fa and (by. Since the dependence on fa and (b ; is only 
through the exponential factors 



e 



±2% g ±2% (16) 
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it is clear that the new subtraction terms integrate to zero after the integration over the helicity 
angles (|), and §j. We therefore modify only the subtraction terms for the real part, but not the 
integrated subtraction terms. We also note that for the method of random polarisations the local 
subtraction terms for the individual helicity eigenstates £+*£^~ and need not be known. 

We would like to add a comment on dimensional regularisation. In the discussion above we 
have assumed for all particles two helicity eigenstates, as if we would work in four dimensions. 
This is justified by the following argument: The subtraction terms are required to match the 
singular behaviour of the integrand in D dimensions. Throughout this paper we use D = 4 — 2e. 
There are several variants of dimensional regularisation. All of them continue the integration 
over the loop momenta and the integration over the phase space of unobserved particles to D 
dimensions. They differ in how they treat the spin degrees of freedom of the particles (and the 
momenta of the observed particles). Conventional dimensional regularisation (CDR) continues 
all spin degrees of freedom to D dimensions, while the 't Hooft-Veltman scheme (HV) continues 
only the spin degrees of freedom of the unobserved particles to D dimensions and keeps the 
ones of the observed particles in 4 dimensions. The four-dimensional helicity scheme (FDH) 
keeps all spin degrees of freedom in 4 dimensions. In the CDR scheme and in the HV scheme 
the unobserved particles may have additional spin degrees of freedom, which are referred to 
as £-helicities. We will discuss the case of a gluon. Let us denote the polarisation vectors 
corresponding to the additional £-helicities by 

(17) 

The spin sum is then given by 

e;X+^X+4^ r ^ 2£) - (18) 

In the four-dimensional sub-space we may replace the sum over the helicities "+" and "— " by 
an integration over the random polarisation E^fy). We thus have 

e^revfo) +4- 2e) *4 _2£) = e+*e+ + e;*£^ +4" 2e) *£^ 2e) + e~ 2 %+*^ + e 2l %*e+. 

The right-hand side is again given as the (D-dimensional) spin sum plus terms proportional to 
£+*£7 and £^*£^~. By construction the subtracted real contribution is integrable in four dimen- 

sions. Therefore we can take the limit D — > 4. In taking this limit all terms proportional to £^ 
disappear. In summary, we have seen that for the extension of the dipole subtraction method 
to the method of random polarisations it is sufficient to consider the four-dimensional helicities 
"+" and "— " only. In the discussion above we have used the gluon polarisations as an example. 
The discussion for the polarisation states of the quarks and the anti-quarks is similar. 

3 Soft and collinear limits 

In this section we discuss the behaviour of \S$\ I 2 in the soft and collinear limits. In the soft 

I n-\-i I 

limit we parametrise the momentum of the soft parton pj as 

Pj = M (19) 
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and consider contributions to |-^+\ | 2 of the order X 1 . Contributions to |-^+\| 2 which are less 
singular than IT 2 are integrable in the soft limit. In the collinear pi\\pj limit we parametrise the 
momenta of two collinear massless final-state partons i and j as 



k 2 

Pi zp - A 



z 2p -n 

Pj = (l-z)p-k ± - f 1 - H . (20) 

1— z2p-n 

Here n is a massless four-vector and the transverse component k±_ satisfies 2pk^ = 2nk^ = 0. 
The four- vectors p, p t and pj are massless: p 2 = pj = pj = 0. The collinear limit occurs for 

k\ — > 0. We consider contributions to |-^+\| 2 of the order |fc^| _2 . Contributions to |^+\| 2 
which are less singular than |fcj_|~ 2 are integrable in the collinear limit. 

The collinear limit, which occurs in massless theories has a generalisation to massive theo- 
ries, called the quasi-collinear limit. In the quasi-collinear limit we consider a splitting ij — > i + j, 
with associated particle masses m;j, m ; and mj. In the quasi-collinear limit we parametrise the 
momenta as 



Pi = zp + k 



_L 



kl+z 2 m 2 j-m 2 n 



z 2p-n 

n \ v fc 2 +(l-z) 2 m 2 .-m 2 n 

Pj = (l-z)p-k ± J - J -- . (21) 

l—Z 2p-n 

Again, n is a massless four-vector and the transverse component k±_ satisfies 2pk^ = 2nk^ = 0. 
The four-vectors p, p, and pj are on-shell: 

P 2 = m fj, P 2 i= m h P 2 j= m2 j- ( 22 ) 

In the quasi-collinear limit we take terms of the order 0(k±), 0(m;y), 0(m ; ) and 0(nij) to be of 
the same order 

0(k ± ) = 0(m ij ) = 0(m i ) = 0(m j ) (23) 

and consider terms which correspond to the order of |&a|~ 2 . Obviously, the collinear limit is a 
special case of the quasi-collinear limit. 

If the emitting particle is in the initial state, the collinear limit is defined as 

Pa = P, 

\ k 2 n 

Pi = {l-x)p + k±- 



1 —x2p-n 
k 2 n 

Pai = xp-k ± . (24) 

x 2p ■ n 
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Here, all particles are massless. In this paper we restrict ourselves to massless incoming partons, 
therefore we do not have to consider the generalisation to the massive quasi-collinear case for 
initial-state partons. 

In the soft limit a Born amplitude with (n + 1) partons behaves as 

lim 4% = gn%(p } )F^ ] . (25) 

Here, g denotes the strong coupling, /j is a scale introduced to keep the coupling dimensionless, 
Pj is the momentum of the soft gluon. The eikonal current is given by 

J" = t T i—- W 

«tl Pi ' Pj 

The sum is over the remaining n hard momenta pi. The colour charge operators for the 
emission of a gluon from a quark, gluon or antiquark in the final state are defined by 

quark : T^A (...qj...) = (T$) A (...qj...) , 

gluon : T^ a SL (...«*...) = (if**) A (...«*...) , 

antiquark : T 9 ^ 9g A (...qj-) = (~T$) & {-qj-) • (27) 

The minus sign for the antiquark has its origin in the fact that for an outgoing antiquark the 
(outgoing) momentum flow is opposite to the flow of the fermion line. The corresponding colour 
charge operators for the emission of a gluon from a quark, gluon or antiquark in the initial state 
are 

quark : {...qj...) = (-7$ A (...qj...) , 

gluon : T^A (../...) = (if ab ) A (.../. 
antiquark : T q ^ qg A (...qj...) = (7$) JZ (...qj..) . (28) 

In the amplitude an incoming quark is denoted as an outgoing antiquark and vice versa. For the 
squares of the colour charge operators one has 

^q-^qg = CF) ^g^gg = Ca- (29) 

We also define the colour charge operator for the emission of a quark-antiquark pair from a gluon 
by 



*(...«*...) = (7g) *(...«*...) (30) 



and 



r Tg^qq = T R- ( 31 ) 
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Ca, Cf and Tr are the usual SU (N c ) colour factors 

N 2 - 1 1 
C A = N C , C F =^-, T R =-. (32) 

In squaring an amplitude we obtain terms proportional to T, • (k ^ i) and terms proportional 
to T?. In order to obtain the dipole structure we would like to re-express T? as a combination of 
terms involving only T, ■ with k ^ i. This can be done using colour conservation. We write 

for ie{q,g} 



■ £ TV TV (33) 

k=\,k^i 



For the splitting g — > qq we write 



= ~ ± ^TrT,. (34) 

k—l,k^i i 

In the following we will often encounter functions which depend on the polarisation of a particle, 
for example in the case of a gluon 

/(e„). (35) 

If the polarisation vector £^ corresponds to a random polarisation with helicity angle § according 
to eq. d8]) then we simply write / (<])) instead of / (e M (<]))), i.e. 

/«>) = f(e M W). (36) 

If on the other hand £„ corresponds to a helicity eigenstate ej; we write 

/(*) = /(«£)■ (3V) 

We use the notation 

/(A) (38) 

to denote either / (<])) or / (X) . This allows for a uniform notation for the polarisations of quarks 
and gluons. 

Squaring the soft limit of the amplitude and using colour conservation one finds 



lim 

pj-X> 



«(0) 

" / Vt-1 



2 ' ' ' ' ■ 

-4na sP 2£ £ £ ^S 0) T,- • (e ; -) (39) 



(=i fe=i 



with 



( P/ • £ ) p, • ej) ( Pi ■ e*) ( Pk ■ £j) + ( Pk ■ £ ) ( Pi ■ ej) 

(Pi-Pj) {Pi-Pj){pi-Pj+Pj-Pk) 
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We will need the (quasi)-collinear limit of the singular soft function Sjj^ (£/) • We find 

lim Sij^Ej) = V ) . (41) 

PiWPi [Pi-Pj) 

The second term of Sy^ (£/) scales like and is therefore not singular enough in the (quasi)- 

collinear limit. Within the method of random polarisations we are in particular interested in the 
terms with a non-trivial dependence on §j. From eq. ([T3l it follows that these terms are given by 



(42) 



Explicitly, we have 
S U,k (0,) 



-2^ 



(43) 



W-ej A-ej [Pre} [Pk-*j +(Pk-*J (Pr*j 



(PrPj)' 



Pi ■ e ; 



Pi ■ el 



(prPj) (Pi-Pj+Pj-Pk) 

Pre/) (Pk-tf) + (pk-ej*) (prtf 
(Pi-Pj) (Pi-Pj + Pj-Pk) 



Let us now turn to the quasi-collinear limit. We discuss immediately the quasi-collinear limit, 
the collinear limit being a special case of the quasi-collinear limit. It will be sufficient to discuss 
the case where all particles are in the final state. The collinear limit in the initial state case can 
be obtained from the collinear limit in the final state by crossing. In the quasi-collinear limit the 
Born amplitude factorises according to 



X T&nll {-,Pi,-,Pj,-) 
Pi\\Pj 



1 



(44) 



where the sum is over all polarisations of the intermediate particle. The splitting amplitudes Split 
are universal, they depend only on the two momenta becoming (quasi)-collinear, and not upon 
the specific amplitude under consideration. The splitting functions Split are given by 

1 _ , 

-u(Pi)$(Pj)u (p), 



Split ^ (p,pupj,huhj) 



Split£ {p,p u pj,hi,hj) 



Splitg^ (p, pi, pj, hj, hj) 



(Pi + Pj) 2 -mfj 
2 U(Pi) -z(pj) Pi- £ x (p)* +e(Pj) • &\p)* Pj-e(Pi) 

^Pi ' Pj L 

-e( Pi ) ■ e\ p )* Pi-e(pj) 



1 u{ P i)4\p)%{ P] ). 



3-Pi-Pj 



(45) 
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Here we used the notation tf^(p) * = e^(p)* y", i.e. complex conjugation is only with respect to 
the polarisation vector. These formulae hold for any polarisation of the particles i, j and (ij). 
We will use these formulae where particles i and j have random polarisations, while the mother 
particle (ij) has a definite polarisation X = + or X = — . There is no point in introducing a random 
polarisation for the mother particle (ij), since eq. (|44]) involves a sum over the intermediate po- 
larisation (or equivalently an integration of the intermediate polarisation). We define the squares 
of the splitting amplitudes by 



[ P {ij)^i+j(PiPhPh h U h j)]^ 



X,X' 



£ ej(p) * Split* * Split*' eJ' (p) for gluons. (46) 

x,x> 



Let us denote by J^f ^ the amplitude, where the polarisation vector of particle (ij) has been 
removed. If particle (ij) is a gluon, \ is a Lorentz index, while in the case where particle (ij) is a 
quark \ corresponds to a Dirac index. With this notation the squared amplitude factorises in the 
collinear limit as 

2 



lim 

PiWpj 



ft. 



(0) 



77+1 



4%a sf ^j4 {0) *T 2 {ij) ^ i+j [P^i+j (p,p h p h h h hj)]„,4' (0) . (47) 



Again we are in particular interested in the terms with a non-trivial dependence on (|), or §j. These 
terms are given by 

P m^i+j (to <k) " I (** M • («) 

Xi ,A, j 

The second term in eq. (1481 corresponds to the splitting functions summed over the polarisations 
of the collinear particles i and j. Subtracting out the sum over the polarisations of two particles 
is an operation which will occur more often in this paper and it is convenient to introduce a 
short-hand notation for that. We define the operator ^ acting on a function f(hi,hj) depending 
on two polarisations by 

!KS(Mj) = f(Mj)-Y,f(b,*.j). (49) 

X.j,A j 

Let us for a moment specialise to the massless case. For the singular part of the splitting functions 
one has then the well-known expressions 



Xi,Xj 
Xi,Xj 



Xi,Xj 



^Pi-Pj 

2 

2Pi-Pj 
+ 
2 



2z 

l-z 



+ (l-e)(l-z) 



2z | 2(1 -z) 



+ 



1 



4(l-£)z(l-z)- 



k 2 



1 



2p r pj 



^ v +4z(l- z )- 



+ 



(50) 
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Note that these equations hold up to terms of order ■ Terms of order 1 / \k± | are integrable in 

the collinear limit. We have defined the splitting functions through eq. (145! ) and eq. (1461) . Through 
these definitions we pick up integrable terms not shown explicitly on the r.h.s of eq. (|50l ). 

The quasi-collinear splittings q—t qg and g — >■ gg have non-vanishing soft limits. We consider 
the case where particle j becomes soft. We find 



(Pi- £ f) (Pi'*]) 

lim [Pq^qg] a Q = -2 Ua(pi)u^pi) 

P.i^° [Pi-Pj) 

(pref) (pi-zj) 
"<-" '" V (Pi-Pj)' 



\imjP g ^ gg ]^ = ^(Pi)*ev{Pi)- (51) 



This shows that the (quasi)-collinear limit of the soft singular function agrees with the soft limit 
of the (quasi)-collinear singular function. To construct the subtraction terms we can therefore 
start from the (quasi)-collinear singular function and supplement it with the terms which are 
singular in the soft limit, but not in the (quasi)-collinear limit. It is convenient to define two 
functions for the terms singular in the soft limit, but not in the (quasi)-collinear limit. We set 



(pi'*j)(pk'Bj) + (pk'ej)(prej) , s / 
l S ^<ig{P>Pi>Pj>Pk,hi,hj)\ a p = , w — — \ MPi)up(pi 

\Pi " Pj) \Pi ' Pj ~r Pj ■ Pk) 



(pi- e j) (Pk ■ ey) + (pk ■ e ■) (pi -ej) 

[ S g^gg\PiPiiPjiPki h i,hj)\ = V7 ■ r E M (pi) 6v(pf) 



{Pi-Pj)(Pi-Pj+Pj-Pk) 

(pj-e*) (p k -£i) + (Pk -e^ipj-Ei) 



%{PjY^{Pj) 



(Pi-Pj) (Pi-Pj+Pi-Pk) 

(52) 



4 The subtraction terms 



In this section we give the additional subtraction terms, which extend the dipole subtraction 
method to random polarisations of the external partons. We recall that the NLO-contribution to 
an observable in the spin-summed case is given by 

(0) NL0 = j (o n+l do R -O ll do A ^+jlo n do v + O n do c + O n Jdo A \. (53) 

n+l n \ 1 / 

We modify this scheme by including an additional subtraction term do A 

(0) NL0 = j O n+l do R - O n (d^ + da^ + J I O n do v + O n do c + O n Jdo A \. (54) 

n+l n \ 1 / 
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The additional subtraction do A term ensures that the expression in the square bracket in the first 
term of eq. (l54l) is locally integrable when used with random polarisations. do A is chosen such 
that 



271 271 



2 / d h / d ^i d<sk = °, (55) 

y J o o 

where (|) ; and (|) ; are the helicity angles of the two unresolved particles. Therefore do A is a 
function, which integrates to zero and can be added to eq. (1531) . Similar to the standard case, do A 
is given as a sum of dipoles: 

^ = EE + E E ^ + E E *? J + E E < 56 ) 

As in the spin-summed case, the dipoles consist of dipole splitting functions sandwiched be- 
tween tree-level amplitudes with n partons. We can use crossing symmetry to obtain all addi- 
tional dipole splitting functions from the final-final case. (The tree-level rc-parton amplitudes are 
evaluated with mapped momenta, and the actual form of this mapping depends on whether the 
partons are in the initial or final state.) In connection with crossing symmetry it is useful to define 
the following operation 

Q : Ei^E*, 

Ui -H- Vj, 

m <H- V(, (57) 

which adjusts the polarisation vector or spinor of the i-th particle from the final to the initial state 
and vice versa. 

4.1 The additional subtraction terms 
4.1.1 Final-state emitter and final-state spectator 

If both the emitter and the spectator are in the final state, the additional subtraction terms are 
given by 

%* = -47ra,A/ 2e (58) 
^ (0) (...,pij,...,p k ,...)* lJ h [Vij !k (pij,pi,pj,p k ,<bi,§j)]£g (0) (...,pij,...,p k ,...) . 

ij 

The functions Vy^ are given for the various splittings by 

Vqigj,k{p,PUPj,PkAAj) = C F H^[P q ^ qg (p,Pi,pj,Mj) +S q ^ qg (p,Pi,Pj,Pk,Mj)] J 
%igj,k PUP jiPh M j) = C A%. [Pg^gg {P, Pi, P j, §U tyj) + Sg->gg (p, PUP j,Pk, <t>i, tyj)], 
V qi q h k(p,PUPj,PkMAj) = T]&.[Pg-+qq{p,PUPj,Mj)]- (59) 
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The operation is defined in eq. (1491) . The mapped momenta pij and pt are defined in the 
massless case by 

/>y =Pi + Pj~ i Pit) = i Pit) J = ; ; ■ (60) 

l -y 1 -y Pi- pj+Pi- Pk+pj- Pk 

In the massive case we use 



X(Q 2 ,m 2 p m 2 ) / Q. pjt \ 2 2 + m2-m 2 



ft = / Pfc~~^2~Q + ^2 2, 

/Me 2 ,(P*+P,) 2 ,m2) v e J 2 Q 



Pij = Q-Ph (61) 
where Q = Pi + Pj + Pk and X is the Kallen function 

X(x,y,z) = x 2 +y 2 +z 2 -2xy-2yz-2zx. (62) 
Eq. (loTI) reduces in the massless limit to eq. (|60t . 

4.1.2 Final-state emitter and initial-state spectator 

If the emitter is in the final state and the spectator in the initial state, the additional subtraction 
terms are given by 

Wf } = -47ra,A/ 2e (63) 

-^ (0) (-,Pij,-,Pa,-)* Wj (PihPUPhPaA^j)]^ *k (0) (-,Pij,-,Pa,-)- 

The dipole splitting function is related by crossing to the final-final case: 

Vfj (Pij, Pi, P j,Pa, <|>i, 4» j) = %a (pij, Pi, P j, ~Pa, <|>i, 4> j) ■ (64) 

The mapped momenta pij and p a are defined by 

Pij = Pi + Pj—0--x)p a , Pa=Xp a - (65) 

The variable x is given by 

Pi-Pa+Pj-Pa-Pi-Pj + j (m 2 j-m 2 -m 2 ) 
x = — , (66) 

Pi-Pa+Pj-Pa 



and reduces in the massless case to 

= Pi-Pa+Pj-Pa-Pi-Pj 
Pi-Pa+Pj-Pa 
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(67) 



4.1.3 Initial-state emitter and final-state spectator 

If the emitter is in the initial state and the spectator in the final state, the additional subtraction 
terms are given by 



(68) 



j4 {0) {...,Paj,..,p k ,...y 



Vk 3 (Paj,Pa,Pj,Pk,§a,§j)\ ^ (-,Paj, -ypk, •••) • 



The dipole splitting function is related by crossing to the final-final case: 

V* ' (Paj,Pa,Pj,Pk,tya,tyj) = C( a ,aj)Vaj,k (~Paj ,~Pa, Ph Pk,$a,$j) ■ (69) 
The operation C is defined in eq. (1571) . The mapped momenta p a j and are defined by 

n v Pk-Pa+Pi-Pa-Pi-Pk ~ m 

Paj=XPa, Pk= Pk+Pi-{l-X)Pa, *= : • (70) 

Pk-Pa+Pi-Pa 

Note that we restrict ourselves to massless initial-state particles. This implies that the masses of 
the particles a, (a j) and j are zero. 

4.1.4 Initial-state emitter and initial-state spectator 

If both the emitter and the spectator are in the initial state, the additional subtraction terms are 
given by 



Lb 



-4%a s p 



2e 



(71) 



i% (0) / ~ ~ \* T £ , j ■ Tfc 



(lh Par-) 



T 2 

aj 



V UJ ' b {paj ,Pa ,Pj ,Pb, <ta Aj ) J *n ^ (p 1 , • • • , Paj , ■ ■ ■ ) • 



The dipole splitting function is related by crossing to the final-final case: 

V ahb (Paj,Pa,Pj,PbAaAj) = C a ,ajV a j,b (~ Paj , ~ Pa, Pj, ~ PbAaA.j) ■ 

In this case the mapped momenta are defined as follows: 

Pa' Pb~ Pi' Pa~ Pi ■ Pb 



Paj=XPa, Pb=Pb, 

and all final state momenta are transformed as 

Pi = Ap/, 
where A is a Lorentz transformation defined by 



Pa-Pb 



(72) 



(73) 



(74) 



A 



(KV + KV)(K V + K V ) 



(K+zy 



+ 2- 



K 2 



K = Pa+Pb~Pj, K = p aj +pb. 



(75) 



Again we consider only the case of massless initial-state particles. Therefore the masses of the 
particles a, b, (a j) and j are zero. 



16 



5 Conclusions 



In this paper we extended the subtraction method for NLO calculations to random polarisations 
of the external particles. We therefore may replace in the computation of the real contribution the 
sum over the helicity amplitudes by a smooth integration over helicity angles. We have derived 
the required additional subtraction terms, such that they match locally the singular behaviour of 
the squared amplitude with random polarisations. These additional subtraction terms integrate 
to zero and modify only the unintegrated subtraction term. They can be used on top of any sub- 
traction scheme, which provides local subtraction terms for the spin-summed amplitude squared. 
Furthermore they are independent of any explicit definition of polarisation vectors or polarisation 
spinors. 
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A Polarisation vectors and polarisation spinors 

We define the light-cone coordinates as 

P+=Po + P3, P-=Po~P3, P±=Pi+ipi, Pi.* = P\-ipi- (76) 

In terms of the light-cone components of a light-like four-vector, the corresponding massless 
spinors (p ± | and \p±) can be chosen as 

\ P+ ) = ( -r ) , \ P -) - ( p 



e 2 e 2 

(P+\ = —7t=t(-P±,P+), {P-\ = —7=={P+iPA-*)i (77) 

v\p+\ V\p+\ 

where the phase § is given by 

p+ = \ p+ \e^, 0<(^<27r. (78) 
If the Cartesian coordinates po, p\, p2 and p^ are real numbers, we have 

\p±)1 = e»{p±\, (p±\ t = e i ^\p±), ^ = ±1. (79) 
Spinor products are denoted as 

(pq) = (p-W+), [qp] = (q+\p-)- (80) 

Let q be a light-like four-vector. We define polarisation vectors for the gluons by 

„+_ (q-\Oy\p-) {q+\o 1! \p+)_ 

% V2(q P ) ' *" V2[pq] ' ^ 
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with = (1,5) and = (1,— a), where d = (01,02,03) are the Pauli matrices. The de- 
pendence on the reference four-vector q drops out in gauge invariant quantities. In numerical 
calculations a convenient choice 11731 for the reference four- vector q for a gluon with momentum 
P/j = (po,Pi,P2,P3) is given by = (po, —p\,—p2, —pi)- With this choice the spinor products 
(qp) and [pq] in the denominator of the polarisation vectors vanish only in the soft limit. In 
the soft limit the numerator vanishes at the same rate. This ensures that the evaluation of the 
polarisation vector is always stable. Under complex conjugation we have 

(«£)*= e;, fe)*=<- (82) 

The reference four- vector q can be used to project any not necessarily light-like four-vector P on 
a light-like four- vector P°: 

The four-vector P b satisfies (P b ) 2 = 0. Let P be a four-vector satisfying P 2 = m 2 . We define the 
spinors associated to massive fermions by 

u± = <F±bi> {If+m) lqTh U± = foiW** 1 {If+m) ' 

v t = tit^ — r(tf-™)kT), v^ = - l ——{q^\(If-m). (84) 

These spinors satisfy the Dirac equations 

(j/-m)u X = 0, u X (p/-m)=0, 

(p ! + m)v X = 0, v X (p/ + m)=0, (85) 



the orthogonality relations 



u X u X = 2m8 U , vV = -2m8 U , (86) 



and the completeness relation 



£«V" = // + m, £vV = //-m. (87) 

X X 

We further have 

u X fu X = 2p^ XX : 

v l fv x = 2p^ xx . (88) 
In the massless limit the definition reduces to 

u ± = v T = (p±\, u ± = v T = \p±). (89) 
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Let us denote the helicity projection operators by 

P+ 



p- = W-y5)=(° n ? )• (90 



2 V V01 
Then 

\p+)(p+\=p+e / , \p-)(p-\=P-rf, (9D 

The combinations \p+) (p—\ and |/?— ) are a little bit more complicated: 

\P+)(P~\ = -j=P+l/lt-(j>,q) = --j=P+lf-(j>,q)l/, 

\P~)(P+\ = -^P-^ + (p,q) = ^P-^ + (p,q)p:. (92) 
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